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1 Introduction 

The gauge hierarchy problem has been one of the main motivations for physics be- 
yond the Standard Model (SM). One of the directions that has emerged, looks at why 
gravity appears weak as compared to the other three interactions of the SM. This 
apparent weakness has been accounted for by the existence of either large extra spa- 
tial dimensions ADD model P or warped extra dimension RS model |2]. In either 
case the fundamental Planck scale could be of the order of a TeV and hence a pos- 
sible explanation of the hierarchy. In both these models only gravity is allowed to 
propagate the extra dimensions while the SM particles are constrained on a 3-brane. 
This consequently leads to Kaluza-Klein (KK) modes in 4-dimensions, which in ADD 
and RS models have distinct KK spectrum and their effective interaction with the SM 
model particles. The experimental signatures of these KK modes have been of intense 
phenomenological activity. With the closure of LEP and the advent of LHC the focus 
now shifts to hadron colliders. 

At hadron colliders, it is important to have a precise knowledge of the parton 
distributions functions (PDFs) to predict production cross sections of both signals and 
backgrounds. These universal PDFs are non-perturbative inputs that are extracted 
from global fits to available data on deep- inelastic scattering (DIS), Drell-Yan (DY) 
and other hadronic processes. They describe the momentum distribution of the partons 
in a proton and various groups have parametrised the PDFs for a wide range of proton 
momentum fraction x carried by the parton and for the center-of-mass energy Q^ at 
which the process takes place. Parametrisation of PDFs to a particular order in QCD 
would involve various theoretical and experimental uncertainties. Recently there has 
been a series of papers [Sj-IOl which for the first time have calculated the next-to- 
leading order (NLO) QCD corrections to various distributions of DY process for both 
ADD and RS model. In [4J, NLO-QCD corrections to dilepton production at hadron 
colliders in the ADD model were presented for the first time, this was extended to 



the RS model in |3]. Further in [Hj we have considered the double differential cross 
section, (Pa /dQ"^ /d cos 9* , for dilepton production in models with large and warped 
extra dimensions. The cos 9* distribution is the one that is actually used by experiments 
and hence is of particular importance. These NLO results would certainly reduce one 
aspect of the theoretical uncertainties as results prior to this calculation were only at 
leading order (LO) in QCD for process involving gravity. 

In j31 El we have quantified the theoretical uncertainties coming from the QCD 
corrections by computing all the processes that enter at NLO level. Unlike the standard 
model contribution to DY, the extra-dimensional models bring in more processes even 
at the LO level. For example, in these models, the gluon initiated process enter at 
the LO level in addition to quark anti-quark initiated process. At LHC, the gluon 
initiated process is more sensitive to factorisation scale compared to quark initiated 
process that necessitated the relevance of NLO computation. It was found in [U E] 
that the NLO corrections are considerably large and the factorisation scale uncertainty 
goes down significantly with these corrections as expected. Our entire analysis is 
model independent because the QCD corrections factor out from the model dependent 
quantities. In jHElin], we used the MRST parton density sets. It is well known that the 
different PDF sets themselves can affect the theoretical predictions and it is important 
to quantify these effects in the observable that could probe new physics. With this in 
mind we have performed a model independent analysis on uncertainties coming from 
the choice of PDF sets in order to make our predictions more reliable. In this paper we 
have looked at the dependence due to various PDFs for the production of dilepton at 
LHC and Tevatron including gravity effects in the ADD and RS models incorporating 
the NLO QCD corrections. The PDF sets used in this study are Alekhin [Zj, CTEQ 
Ej and MRST jHj . Dependence on PDF sets is also compared with experimental errors 
that enter the parametrisation of the PDF, which are now available to NLO QCD 
[inillll. For this purpose we used the MRST distribution [Tl] as a typical case. The 



dependence of factorisation scale fip and renormalisation scale fiR in going from LO 
to NLO is also studied. 

For the ADD model we also study the dependence on UV cutoff of the KK mode 
sum by keeping the UV cutoff different from the scale of the model Ms- This prescrip- 
tion allows us to study the cutoff dependence. The dependence of the cross section on 
the number of extra dimension as a result of this prescription is similar to that of the 
real graviton production case. The ADD model is a low energy effective theory valid 
below the scale Ms, it is conventional to equate the cutoff to the scale of the effective 
theory ^^. 

The plan of the paper is as follows. In section 2 we describe the extra dimension 
models studied and describe the parameters of the models. In section 3 we discuss the 
theoretical uncertainties viz. the PDF uncertainties, renormalisation and factorisation 
scale dependences. The improvement of the scale dependence in going from LO to 
NLO is also discussed. In section 4 we look at the dependence due to the error on the 
data and estimate the experimental error for a few observable using the MRST PDF. 
Finally we summarise our results in section 5. 

2 Extra Dimension Models 

Extra dimension models that allow gravity to propagate the extra dimensions would 
in 4-dimensions have KK modes which couple to SM particles through the energy 
momentum tensor. The Feynman rules of the KK mode interactions with the SM 
fields are given in J21 ^] ■ Due to the different methods of compactification of the extra 
dimensions in ADD and RS models, their KK spectrum are very distinct. Experimental 
signature of extra dimensions would correspond to deviation from SM predictions due 
to the virtual exchange of KK modes or direct production of KK modes at a collider. 
In the ADD case, there is a tower of KK modes which are almost degenerate in 
energy and a sum over these KK modes gives an observable effect. In the case of 



dilepon production, in addition to the SM photon and Z production modes, one has to 
take into account virtual KK modes. Performing the sum over the virtual KK modes 
leads to an integral which has to be regulated by an UV cutoff. The propagator after 
the KK mode summation becomes 



where n = y/lGn/Mp is the strength of the gravitational coupling to the SM particles, 
rrin the mass of KK modes, d is the number of extra dimensions and Ms is the scale of 
the A + d dimensional theory. The summation over the non-resonant KK modes yields 
I{Ac/Q) jni- Conventionally the UV cutoff A^ is identified with the scale of the extra 
dimension theory Ms, which simplifies the expression giving a mild dependence on the 
number of extra dimensions ^21 ^j • 

In this analysis, we have kept the cutoff Ac different from Ms ^. Note that the 
summation of KK modes in Eq. (Q), modifies the Mp suppression to Ms suppression. 
The ADD model is a effective low energy theory valid below the scale Ms, for con- 
sistency it is essential to satisfy the condition Q < Ac < Ms- The parameters of 
the ADD model are Ms the scale of the 4 + d dimensional theory and d the number 
of extra spatial dimensions. If Ac 7^ Ms then there is an additional parameter. We 
have studied the dependence of the cross section on the cutoff Ac = aMs and varied 
a = 0.7 — 1. In Fig. [T^ we see that the cross section decreases as we lower the cutoff 
Ac. The corresponding K-factor also decreases for lower cutoff Fig. ^p. Dependence of 
the cross section on the number of extra dimensions d is shown in Fig.^ for Ac = Ms, 
the cross section decreases as d increases. Reducing Ac decreases the cross sections 
and if d is increased it brings down the cross section much faster Fig. ^. Large extra 
dimension searches in the dimuon channel at the Tevatron [15^ have put bounds on 

Ms in the range 0.8 - 1.27 TeV. For the analysis, they have used the double differential 

§ Effects of the various UV cutoff methods on the low scale quantum gravity model have been 
discussed in |14| . 
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cross section with respect to invariant mass and the cos 9* ^H] • 

In the RS model, the gravity propagate one extra dimension which is warped by 
an exponential factor exp(— vr/cL), where L is the compactification length and k is 
the curvature of the AdS^ space-time. The parameters of the RS model are itiq = 
kexp{—TTkL) which sets the mass scale of the KK modes and Cq = k/Mp the effective 
coupling. The higher KK modes have enhanced coupling to SM particles due to the 
warp factor and decouple from the zero mode, which is as usual Mp suppressed. RS 
KK spectrum is distinct from the ADD case and hence the summation of the KK 
modes that contribute to the virtual process would also be different. The function 
P(Q^) in the KK mode propagator results from summing over the resonant KK modes 
and is given by 

^(Q') = T.7rr-Tirr-wir^-2^ (2) 



- Q2 _ M2 + tMnVr 
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where M„ are the masses of the individual resonances and F^ are the corresponding 
decay widths. The graviton widths are obtained by calculating their decays into final 
states involving SM particles. A is defined as 
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where Xg = Q/rriQ. We have to sum over all the resonances to get the value of X{xs), 
which is done numerically for a given value of Xg- Searches for the RS KK modes at 
Tevatron in the dielectron, dimuon and digamma channel ^\ have yielded a lower 
hmit between 250 - 785 GeV depending on the coupling to the SM particles. 

3 Theoretical uncertainties 

In the QCD improved parton model the hadronic cross section can be expressed in 
terms of pertubatively calculable partonic cross sections denoted by a"^(r, (5^,/iF) 
convoluted with appropriate non perturbative partonic flux $a;,(r, /i^) at a factori- 
sation scale fip. The subprocess cross section is a perturbative expansion in the strong 



coupling constant asifJ^R) and is calculated order by order in as- Here hr is the renor- 
malisation scale and r = Q'^/S is the DY scaling variable. In perturbative QCD, the 
unknown higher order corrections and the scale uncertainties are strongly correlated. 
The factorisation of mass singularities from the perturbatively calculable partonic cross 
sections leads to the introduction of factorisation scale hf in both non-perturbative 
partonic flux ^abif^p) as well as the finite partonic cross sections daab{x, fip)- Even 
though the choice of the scale is guided by the hard scale of the problem, the exact 
value does not come from the theory. The PDFs and partonic cross sections satisfy 
renormalisation group equations such that the hadronic cross section is independent 
of the factorisation scale fip- In addition to the factorisation scale, the partonic cross 
sections are dependent on the renormalisation scale hr. The choice of the scale is again 
arbitrary. Even though this is an advantage to choose appropriately to do perturbative 
calculations, it also introduces theoretical uncertainties through the size of unknown 
higher order corrections. Usually, one chooses this scale such that the perturbative 
methods can be applied and then computes higher order corrections sufficiently such 
that the exact choice of this scale becomes almost immaterial. Gravity couples to the 
SM fields via its energy momentum tensor, and the calculations are done in the high 
energy limits where masses of the SM particles are ignored. Only parameter that re- 
quires UV renormalisation is the strong coupling constant, because of this we have the 
following expansion for the mass factorised partonic cross section: 

oo 

daab{z, ^D = Y, al{fil)da^^{z, /i|, fi^) , (4) 

where the coupling constant satisfies standard renormalisation group equation. Since 
we are only interested in the NLO order corrections, the Altarelli-Parisi kernels P^'^^z), 
P^^\z) and the coefficients /3o, Pi are sufficient for our analysis. The scale uncertainties 
come about from the truncation of the perturbative series. Unlike the perturbatively 
calculable partonic cross sections, the PDFs being non-perturbative in nature are ex- 
tracted from various experiments. These are fitted at a scale of the experiments and 
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then evolved according to the AP evolution equations to any other relevant scale. They 
are not only sensitive to experimental errors but also to theoretical uncertainties that 
enter through the partonic cross section calculations and the splitting functions that 
are known only to certain orders in strong coupling constant in perturbative QCD. 
Here, we mainly concentrate on the uncertainties coming from PDFs in detail and 
quantify their impact on the new physics searches in extra dimensional models. 

3.1 PDF uncertainty 

We first focus on the uncertainties coming from different PDF sets. The parton flux 
factor for both LHC and Tevatron would give an idea as to which component would 
be dominant in the kinematical region of interest. This flux factor enters the cross 
section. The gluon flux is clearly much larger in the kinematical region of interest at 
LHC and for Tevatron the qq flux is the dominant contribution. 

In the context of extra dimension theories we consider the dilepton production at 
LHC and Tevatron for both large extra dimension and warped extra dimension models. 
The process of interest is -Pi(pi) + P2{P2) -^ f^^ih) + fJ'^ih) + X{Px), where Pi and 
P2 are the incoming hadrons, /i^ are the flnal lepton pair and X the flnal inclusive 
hadronic state. The dilepton in these models could also be produced from the exchange 
of a KK mode in addition to the usual SM gauge boson exchange. Hence at LO itself 
the gg subprocess could contribute to the dilepton production via a KK mode exchange 
in addition to the qq subprocess. 

For both new physics searches and precision SM physics at hadron colliders it is 
essential to understand the uncertainties associated with PDFs. We essentially study 
to what extent the cross sections depend on the various PDFs viz. Alekhin [2j , CTEQ 
[H] and MRST jH]. In the Table 1, we have tabulated the particular PDF that is chosen 
for the study and also the corresponding Aqcd parameter that is used to determine 



the strong coupling a^. ^ 



LO 


NLO 


PDF 


Aqcd (GeV) 


PDF 


Aqcd (GeV) 


MRST2001 LO 
CTEQ6L 


0.220 
0.326 


MRST2001 NLO 

CTEQ6M 


0.323 
0.326 



Table 1: The PDF set used in the analysis along with the respective Aqcd- 



These groups perform a global analysis of a wide range of DIS and other scattering 
data to get best fits to a particular order in QCD. Though all these parametrisation 
satisfies the general constraints, they could differ from each other. This is expected as 
PDFs are not by itself physical quantities and are extracted subject to experimental 
and theoretical uncertainties and various assumptions and initial conditions used by the 
different groups. Differences among the various PDFs would translate as uncertainties 
on the physical observable. 

To NLO in QCD for various PDFs, we now present the comparison plots for the 
following differential distributions 

(fa d?a 



da 
dQ' 



(5) 



dQ dY ' dQ dcos9* 

We would look at the invariant mass distribution Q, the double differential cross section 

with respect to Q and rapidity Y and the double differential cross section with respect 

to Q and cos 6*. The angle 6* is the angle between the final state lepton momenta 

and the initial state hadron momenta in the c.o.m frame of the lepton pair. The 

corresponding K-factor which is the ratio of NLO to LO of the above distributions 

are also plotted for the various PDFs. For the double differential cross section we fix 

the invariant mass Q in the region of interest of extra dimensions and plot the cross 

section with respect to rapidity Y and cos^*. The first two distributions in Eq. (0) are 

cos 6* integrated distributions and hence are independent of the interference between 
^In the case of Alekhin the PDF itself generates the value of as and is hence not tabulated. 



the SM background and the low scale gravity effects 4]. The double differential with 
respect to Q and cos 9* would contain the interference terms, but numerically it is not 
very significant [6 . Consequently even for the cos^* distributions we can express the 
K factor of the model involving both SM and gravity as 

where K^^ is the K factor of the pure gravity part. We have introduced a quantity 
K'^^^ , defined as the ratio of the LO distribution of gravity to SM, given by 



K^'\Q) 



dait^Q) 



1 -1 



dQ 



(7) 



dQ 

The behaviour of K^^\Q) is governed by the competing coupling constants of SM and 
gravity and the parton fluxes involved. Basically the factor K^^^ is an indicator as to 
the source of the total K^^^^^^^^-iactor. K^^\Q) as a function of Q rises much faster 
for LHC than Tevatron and reaches 1 much earlier. Since the gg subprocess contributes 
at LO itself for the gravity mediated process, the gravity effects are much larger at the 
LHC where the gluon flux is much larger. This would also result in larger K-factor for 
the process at LHC at large Q where the gravity contribution dominates. At Tevatron 
since the gluon flux is smaller the K-factor is similar to the SM K-factor. 

For both ADD and RS models the signal for new physics is the excess of events 
in the total cross section or the various distribution over the SM background. If we 
restrict ourself to these extra dimensional models, the signal is due to the effect of the 
KK modes and can not be mimicked by the SM. We would like to emphasise that we 
are not analysing the existing Tevatron data to extract bounds on the ADD and RS 
parameters, which would need a full hadron-level simulation, but estimate the various 
uncertainties to NLO in QCD by choosing typical representative values for the ADD 
and RS parameters. 

We begin with the ADD model wherein we have chosen d = "i and Ms = 2 TeV. In 
Fig. 12^^ the cross section is plotted as a function of the invariant mass Q of the dilepton 



at LHC for the various PDFs. There is only a mild dependence on the difference in 
the PDFs, but when plotted for the corresponding K-factor then the PDF dependence 
is larger for both low and high values of Q, Fig. ^p. At low Q it is the SM part which 
contributes to the K-factor while at high Q it is the beyond SM effects that contribute. 
At low Q where the K-factor is due to SM part, MRST and CTEQ are similar, while 
Alekhin is smaller. At large Q the K-factor is due to the gravity part and here CTEQ 
is larger. 

For the double differential cross section with respect to invariant mass distribution 
and rapidity Y Fig. EK, we have plotted as a function of rapidity Y for a fixed Q = 0.7 
TeV. Only in the central rapidity region do the PDFs differ, with MRST being the 
dominant while CTEQ is the smallest. The K-factor is quite large at the central 
rapidity region and would range from 1.5 - 1.6 depending on the PDF used. The general 
behaviour of the K-factor is similar for MRST and Alekhin. At large rapidities y = ±2 
the K-factors are quite different with Alekhin being 1.25 while CTEQ the largest is 
1.45. For Q = 0.7 TeV the K factor is large which we can see from Fig. |2b, wherein 
the dominant contribution is from the gravity mediated gg initiated subprocess. 

In Fig. Efc we have plotted the double differential cross section with respect to Q 
and cos 6** as a function of cos 6'* for a fixed Q = 0.7 TeV. MRST gives the largest and 
CTEQ the least with Alekhin being a central value in the spread. The difference exists 
for the full range of cos 6'*. The SM background has a different cos 6** dependence. The 
interference of the SM and the gravity effect is not zero for the cos 6* distribution but 
does not contribute significantly. The K-factor for central cos 6'* = region is about 
1.52 but differ with PDFs as cos^* -^ ±1, Fig. Eli. Since there is no gg initiated 
process in the SM background to NLO the K-factor is much smaller. 

In the RS model we have chosen the mass of the first KK mode Mi = 1.5 TeV and 
the coupling cq = 0.01. In Fig. |3^ we have plotted invariant mass distribution of the 
dilepton in the RS model. At the KK mode resonances the cross section differs from 
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the SM cross section, but the dependence on the PDFs are very mild. In Fig. EJd the 
corresponding K-factors are plotted for various PDFs. The behaviour of the K-factor 
of the RS model can be understood with the help of Eq. (|6I7|) . It is only in the RS 
graviton resonances region that K^^^ is large and hence the i^-factor is dominated by 
the K'^^ factor. In the off resonance regions it is the K^^^ which contributes. There 
is a wide difference in the K-factor more in the second peak and even off peak where 
the effect is mainly SM. This is due to the high Q value that is chosen in the RS case. 

For the double differential with respect to rapidity and invariant mass, in Fig. Et 
we have plotted it for the rapidity range of LHC for Q = 1.5 TeV, which is the region 
of the first RS KK mode. It is only in the resonance region that the effects of RS 
are visible. Here there seems to be a clustering of PDFs but for CETQ in the central 
rapidity region. In the central rapidity region the K-factor varies from 1.6 - 1.75 Fig.Eb- 
In the first RS KK resonance region at Q = 1.5 TeV the gravity dominates and hence 
the K-factor is large (Eq. ©). Beyond the central rapidity region Y = the K-factor 
dependence on PDFs is substantial. In Fig. Et we have plotted the double differential 
with respect to cos 6* for Q fixed at the first resonance. The cross section is largest for 
cos 6* = and MRST is the largest among the PDFs. The K-factor in Fig. Eli is about 
1.65 for wide range of cos 6'* for Alekhin and MRST but for CTEQ it varies between 
1.7- 1.8. 

In the above we discussed the extra dimension effects at LHC, now we look at the 
Tevatron. For the ADD case, in Fig. IHk we have plotted the invariant mass distribution 
for the various PDFs. The spread due to the various PDFs over the Q^ range is not 
too large. Only at large Q there is some deviation from the SM result which is plotted 
in Fig. int- The K-factor for the Q distribution for the various PDFs are plotted in 
Fig. inb, which are in tune with the SM K-factor at the Tevatron. In Fig. IHfc we have 
plotted the PDF comparison plot for the rapidity distribution at Q = 0.7 TeV. CTEQ 
and MRST plots are very similar while Alekhin is larger in the central rapidity region. 
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In the y = region, the K-factor for CTEQ is about 1.1 while for MRST and Alekhin 
it is about 1.2, which is in the range of the SM K-factor, Fig. IHll. 

For the RS model the PDF comparison plots are given in Fig. [7| In Fig. [7^ we have 
the invariant mass distribution and the deviation from the SM is only in the resonance 
region. The PDF dependence is very mild. In the first resonance region the K-factor 
(Fig. IZb) is dominated by K'^^ at Q = 0.7 TeV but at Tevatron this value is not too 
different from the SM K-factor. In Fig. ^ the cos 9* distribution at the first resonance 
region is plotted, CTEQ and MRST overlap while Alekhin is larger over a wide range 
of cos^*. The K-factor Fig. [2|i is in the range of the SM K-factor. 

3.2 Renormalisation/Factorisation scale uncertainties 

In Fig. IHK we have plotted the double differential (Pa/dQdY in the Y range for LHC 
energies for a fixed Q = 0.7 TeV. The dependence of cross section on ^r comes from the 
strong coupling constant at NLO and so at LO there is no iir dependence. At NLO /xr 
dependence for the Y distribution is plotted for the iir range 0.5 Q < fiR < 1.5 Q. The 
fiR spread is largest in the central rapidity region and would only reduce at the NNLO 
order level when the /ir dependences would be compensated for by the dependence 
coming from the coefficient functions. In Fig. |Hb we have plotted the K-factor for SM 
and SM+GR and see how it dependence on fiR. The uncertainties due to fiR is much 
larger when the gravity is included. The percentage spread is of the order of 3.5 % 
which is comparable to the fip spread at NLO. 

In Fig.iniwe have plotted Y distribution and its K-factor for ADD and RS model at 
a fixed Q = fiR. The ftp variation is studied by varying ftp in the range 0.5 Q < fip < 
1.5 Q. We see that for both the ADD and RS model in going from LO to NLO in QCD, 
the uncertainties due to ^p variation considerably get reduced. The spread of K-factor 
with fip is much smaller for the SM as compared to SM-I-GR. This certainly indicates 
need to go beyond NLO. Similar trends are observed for the cos 6'* distribution plotted 
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Distributions 


Tevatron 


LHC 


LO 


NLO 


LO 


NLO 


ADD 


cPa/dQdY 
<fa/dQdcos9 


22.8 
24.2 


7.4 
8.2 


9.5 
10.9 


3.5 
3.8 


RS 


d^a/dQdY 
d'^a/dQdcosO 


23.2 

24.2 


7.7 
8.0 


18.7 
18.4 


6.9 

6.8 



Table 2: Percentage spread as a result of factorisation scale variation in the range 
0.5(5 < /^F < 1.5Q. For the ADD case Q = 0.7 TeV. For the RS first resonance region 
Q = 1.5 TeV for LHC and Q = 0.7 TeV for Tevatron. 

in Fig. Uni 

In Table 12] we tabulate the percentage spread of the factorisation scale fip depen- 
dence in the range 0.5Q < ^f ^ 1.5Q for the LHC and Tevatron. On the average at 
LHC and Tevatron, the percentage spread of the scale variation get reduced by about 
2.75 times in going from LO to NLO. 

4 Experimental Uncertainties 

In addition to the theoretical uncertainties that we have described in the previous 
section, there are uncertainties due to errors on the data. Various groups have studied 
the experimental errors and have estimates of the uncertainties on the PDFs within 
NLO QCD framework ^0] E]- Now that NLO QCD results are also available for 
extra dimension searches ^ for the dilepton production, we consider some of the 
distributions and estimate the uncertainties due to the experimental error. In Fig. ITTk 
we have plotted the error band for the MRST 2001 PDF ^T] in the ADD model for 
the dilepton invariant mass distribution at LHC. This error band is comparable to the 
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spread associated with the different set of PDFs as given in Fig. ^.. At Q = 1 TeV 
the percentage of experimental error is 7.5 % for SM + GR while the pure SM error 
is about 3.3 %. For the RS case at LHC in the first resonance region at Q = 1.5 TeV 
the experimental error is about 12.8 %. At Tevatron the ADD model experimental 
error is 7.4 % at Q = 1 TeV. The experimental error for this distribution for the 
central rapidity region is about 3.5 % and is indicated in the Fig. Illb. In general the 
experimental error increases with the increase in Q. 

5 Conclusions 

We have studied the impact of various parton density sets at next to leading order in 
strong coupling constant a^ in QCD to one of the most important processes, namely 
Drell-Yan production of dileptons in hadron colliders such as LHC and Tevatron. This 
process can probe the physics beyond SM through exchange of new particles that these 
theories predict. At hadron colliders, the precise measurement of DY production cross 
sections is possible. In this context, we have studied the theories of extra dimensions 
such as ADD and RS which attempt to explain gauge hierarchy problem in SM. We 
have discussed various theoretical uncertainties that enter through renormalisation, 
factorisation scales and the parton density sets. We have quantified the uncertainties 
coming from various parton density sets using the recent results on NLO QCD correc- 
tions to parton level cross sections and recent PDF sets that take into account various 
theoretical and experimental errors. Our entire analysis is model independent thanks 
to the factorisation of QCD radiative corrections from the model dependent contribu- 
tions. More precisely, our findings are independent of the finer details of the model as 
they factor out from the rest. We find that the K- factor for various observable depends 
on the choice of PDFs. 
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Figure 1: (a) Invariant mass distribution is plotted for various values of the cutoff 
Ac = aMg in the ADD model, (b) The corresponding K-factor. (c) Invariant mass 
distribution as a function of the number of extra spacial dimension d for A^ = Ms TeV 
at LHC (d) The same plot as (c) for Ac = O.TM^. 
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Figure 4: (a) The invariant mass distribution of dilepton pair production at LHC in 
the RS model for various PDfs. (b) The corresponding K-factor for the various PDFs. 
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Figure 5: (a) The double differential cross section with respect invariant mass and 
rapidity for various PDFs in the RS model at Q = 1-5 TeV, the region of first resonance, 
(b) The corresponding K- factor as function of rapidity at Q = 1.5 TeV. (c) In the 
region of first RS resonance, the double differential with respect to invariant mass and 
angular distribution of the lepton is plotted for the various PDFs at LHC (d) The 
corresponding K-factor for the various PDFs. 
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Figure 6: ADD model at Tevatron for the various PDF set, we plot in (a) the invariant 
mass distribution. In (b) the corresponding K-factor. (c) The double differential with 
respect to Q and Y is plotted for a fixed Q = 0.7 TeV and for the Y range of Tevatron. 
In (d) the corresponding K factor is plotted. 
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Figure 7: RS model at Tevatron for the various PDF set, we plot in (a) the invariant 
mass distribution. In (b) the corresponding K-factor. (c) The double differential with 
respect to Q and cos6'* is plotted for a fixed Q = 0.7 TeV and for the cos 6'*. In (d) 
the corresponding K factor is plotted. 
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Figure 9: (a) Factorisation scale dependence for the double differential cross section 
as a function of rapidity for LO and NLO for factorisation scale in the range fip = 
0.5 Q - 1.5 Q. (b) SM and SM+GR K factor for ADD rapidity distribution in the 
same variation of fip. In (c) the RS distribution at Q = f^R = 1.5 TeV in the region of 
first resonance, (d) The SM and SM+GR K factor for RS rapidity distribution. 
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Figure 10: (a) Factorisation scale dependence for the double differential cross section 
as a function of cos^* for LO and NLO for factorisation scale in the range [ip = 
0.5 Q - 1.5 Q. In (b) we have plotted the SM and SM+GR K factor for ADD at 
Q = I^R = 0.7 TeV. In (c) the RS cos^* distribution for LO and NLO in the same 
range of iip- (d) The SM and SM+GR K-factor dX Q = iir = L5 TeV, the region of 
first resonance. 
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Figure 11: The experimental error on the MRST PDF at LHC in the ADD model 
for (a) The invariant mass distribution and (b) The rapidity distribution for a fixed 
g = 0.7 TeV. 
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